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Abstract 

Sufficient conditions for wave breaking are found for the short-pulse equation de- 
scribing wave packets of few cycles on the ultra-short pulse scale. The analysis relies on 
the method of characteristics and conserved quantities of the short-pulse equation and 
holds both on an infinite line and in a periodic domain. Numerical illustrations of the 
finite-time wave breaking are given in a periodic domain. 

1 Introduction 

The short-pulse equation, 

(1.1) utx = u + l{u^)xx, t>0,xe'R, 

b 

is a useful and simple approximation of nonlinear wave packets in dispersive media in the 
limit of few cycles on the ultra-short pulse scale [U [TT]. This equation is a dispersive 
generalization of the following advection equation 

(1.2) ut = ^u^Ux, t > 0, X G M. 

According to the method of characteristics, the advection equation (jl.2p exhibits wave 
breaking in a finite time for any initial data u{x,0) = uo{x) on an infinite line if uo{x) is 
continuously differentiable and there is a point xq G M such that uq{xq)uq{xo) > 0. This 
follows from the implicit solution 

uic,t) = uoio, x(e,t) = e-^n2(e)t, t>o, ^EM, 

for any given uq{x) £ C^(R.). We say that the finite-time wave breaking occurs if there 
exists a finite time T £ (0, oo) such that 

(1.3) limsupu{x, t)ux{x, t) = oo, while limsup |ii(3;, t)| < oo. 
For the simple advection equation ()1.2p . 

T = inf I \ : no(e)n'o(0 > o| . 
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In view of this result, we address the question if the dispersion term d~^u in the short-pulse 
equation (II. ip can stabilize global dynamics of the advection equation ()1.2p at least for 
small initial data. Local well-posedness of the short-pulse equation on an infinite line was 
proven in [llj . 

Theorem 1 (Schafer & Wayne, 2004). Let uq G if^(M). There exists aT > such that 
the short-pulse equation admits a unique solution 

u{t) G C{[0,T),H^{R))nC\[0,T),H\R)) 

satisfying n(0) = uq. Furthermore, the solution u{t) depends continuously on uq. 

To extend local solutions into a global solution, Pelinovsky & Sakovich [9] used the 
following conserved quantities of the short-pulse equation: 

(1.4) £^0 := / u^dx, 

(1.5) := [ (yrT^-l)dx= / ^ -dx, 

(1.6) E2 := [ VTT^, 

JR 

If u{t) is a local solution in Theorem [H then Eq, Ei, and E2 are bounded and constant in 
time for all t G [0,T). Global well-posedness of the short-pulse equation on an infinite line 
was proven in [91 for small initial data in satisfying 

(1.7) 2Ei + E2< Wu'qWI^ + HWh < 1. 

The condition p.7p can be sharpen using the scaling transformation for the short-pulse 
equation (jl.ip . Let a G M+ be an arbitrary parameter. If u{x,t) is a solution of (jl.ip . then 
u{x, t) is also a solution of (jl.ip with 

(1.8) X = ax, i = a^^t, u{x,i) = au{x,t). 
The conserved quantities transform as follows: 

El = f Y^l+^- 1) = « / + -Adx = aEi, 

Jr (1 + n|)5/2 (1 + n2)5/2 

Finding the minimum of 2Ei + E2 = 2aEi + a-^E2 in a gives a sharper sufficient condition 
on global well-posedness [9]. 

Theorem 2 (Pehnovsky & Sakovich, 2009). Let uq G H'^{R) and 2y/2E^ < 1. Then the 
short-pulse equation II 1-1]) admits a unique global solution u{t) G C(R-|-, i/2(M)) satisfying 
u{0) = no- 



5t 



dx 



(1 + ^2)5/2 



dx. 
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Theorem [2] does not exclude wave breaking in a finite time for large initial data and 
this paper gives a proof that the wave breaking may occur in the short-pulse equation 
(jl.ip . Negating the sufficient condition for global well-posedness in Theorem [2l a necessary 
condition for the wave breaking follows: the wave breaking may occur in the short-pulse 
equation (jl.ip with the initial data uq £ only if 2\j2EiE2 > 1. We shall find a sufficient 
condition for the wave breaking in the short-pulse equation (jl.ip . 

Unlike the previous work in [9], we will not be using conserved quantity E2 but will rely 
on the conservation of Eq, Ei, and the energy 



dx. 



(1.9) E.,:= [ [d-\)'-l-r 

Here d~^u is defined from a local solution u by 

u{x',t)dx' = — u{x',t)dx' = \ — I u{x',t)dx', 

00 Jx ^ \J ~oo Jx J 

thanks to the zero-mass constraint Jj^u{x,t)dx = for all t G (0, T). (Note that the initial 
data no does not have generally to satisfy the zero-mass constraint J^UQ(x)dx = 0.) 

Thanks to the Sobolev inequality, ||u||^4 < C||n||j|/i for some C > 0, the quantity E^i is 
bounded if u G H'^{M) D i^^^M), where H"^ is defined by its norm 

ll^^lliz-i := \\dx^u\\L2- 

(Note that if u G H^{R)r\H-^{R), then J^u{x)dx = is satisfied.) Our main result on the 
wave breaking on an infinite line is formulated as follows. 

Theorem 3. Let uq G if^(M) H H^^{R) and T be the maximal existence time of Theorem 
m Let 

F, := ^{Ef + {8EoE,+Ef)'/'y^\ 

1 / 1 2^^/2 

Fq :— —7= Eq + £'_! -I- —EoFi 



V2 V 12 
and assume that there exists xq G M such that uo{xo)u'q{xo) > and 



/ ^2 X 1/3 / ^ 1/2 

either |u'o(xo)| > f ^ j , |uo(xo)| |n'o(xo)|' > Fi + f 2Fo|n'o(xo)|=^ - -Ff 

or 1^0*^^°)' - (4^) ' ho(a;o)||no(xo)P > Fi. 



Then T < 00, so that the solution u{t) G C{[0,T), H'^{R) n H-^{M)) of the short-pulse 
equation (CUP blows up in the sense of 

limsupu(a;, t)'U2:(x, t) = 00, while lim\\u{-,t)\\]^oo < Fi. 
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Remark 1. The quantities Fq and Fi in Theorem\^ can he defined by 
Fi := sup \\u{-,t)\\Loo, Fq := sup \\d~^u{-,t)\\L-<' ■ 

te[0,T) tG[0,T) 

Note that the scaling transformation il.8\) gives 

\\u{;i)\\L^ = a\\u{;t)\\L^, \\d^^u{;i)\\L^ = a^\\d-\{;t)\\L^ , 

so that the sufficient condition of Theorem with new definitions of Fi and Fq is invariant 
in a. We note, however, that, while the bound on Fi in Theorem scales correctly as 
Fi = aFi, the bound on Fq is not correctly scaled in a because 

This is an artefact of using Sobolev embedding in Lemma\^ below. 

To prove Theorem[3l we shall adopt the method of characteristics and proceed with apriori 
differential estimates. Our results remain valid in a periodic domain, where Theorem [3] is 
replaced by Theorem S] below. The technique of characteristics and apriori differential 
estimates has been applied for wave breaking in other nonlinear wave equations, see [21 [3l 
HI El m [71 [12] for an incomplete list of references. 

We emphasize that, unlike the previous work in [9], we avoid using a transformation 
between the short-pulse equation (II. ip and the integrable sine-Gordon equation in char- 
acteristic coordinates. Our proof of the wave breaking for the short-pulse equation (jl.ip 
does not suggest that there exists a similar wave breaking for the sine-Gordon equation, it 
is rather the breaking of the coordinate transformation between the two equations. On a 
similar note, we do not use the integrability properties of the short-pulse equation (II. ip such 
as the Lax pair, the inverse scattering transform method, the bi-Hamiltonian formulation, 
and the existence of exact soliton solutions. 

The article is constructed as follows. The proof of Theorem [3] is given in Section 2. 
Section 3 reports extension of Theorem[3]to a periodic domain. Section 4 contains numerical 
evidences of the finite-time wave breaking in a periodic domain. 



2 Wave breaking on an infinite line 

Let us rewrite the Cauchy problem for the short-pulse equation on an infinite line in the 
form 



(2.1) 



ut = "^wux + d^^^u, X G M, t > 
m(x,0) = no(x), e M, 



where d~^u := f^^u{x',t)dx'. In what follows, we use both notations u{t) and u{x,t) 
for the same solution of the Cauchy problem (|2.ip . Local existence of solutions with the 
conservation of E-i and Eq is described by the following result. 
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Lemma 1. Let uq £ i?^(M) n H ^(M), s > 2. There exist a maximal time T = T{uq) > 
and a unique solution u{x, t) to the Cauchy problem 12. 1\) such that 

u{t) e C{[o,T),H-%R)nH-^{w))nC\[o,T),H'-\R)) 

satisfying u{0) = uq. Moreover, the solution u{t) depends continuously on the initial data 
Uq and the values of E^i, Eq, and Ei in {1-4^ , and \1.9^) are constant on [0, T). 

Proof. If no G ^''(M) n H-'^{R), s > 2, then d'^UQ £ H'+'^{R), so that J^uo{x)dx = 0. By 
the theorem of Schafer & Wayne there exists a solution 

u{t) e C{[0,T),H'{R))nC\[0,T),H'-^' 

of the short-pulse equation (jl.ip . so that 

d-\{t) := ut - \u^u, e C{{Q,T),H'-^ 

Therefore, u{t) G C([0, T), i?^(M) n iJ"i(]R)) in view of continuity of ||ii||^-i at t = 0. 
Because / E iJ'*(R), s > 1 implies lim|^|^oQ f{x) = 0, the zero-mass constraint holds in the 
form 



(2.2) 

Let us define 



/ u{x,t)dx = 0, t£ [0,T). 



By the zero-mass constraint (I2.2p and uniqueness of the solution u{t) for any t G [0, T), we 
obtain 

lim d~^u{x,t) = 0, t G [0,r). 

|a;|— >oo 

Using balance equations for the densities of E^i, Eq, and Ei, we write 



1 

12 



-u 



1 

36' 



-u 



l + ui 
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Integrating the balance equation in a; G M for any t G [0, T), we complete the proof that 
-E-i, Eq, and Ei are bounded and constant on [0, T). □ 

Remark 2. The maximal existence time T > in Lemma{l\is independent of s >2 in the 
following sense. If uq G H^{R) n H^' (R) n for s,s' >2 and s / s' , then 

u{t) G C{[0,T),H'{R)nH-\R))nC\[0,T),H'-^{ 

and 

u{t) G C{[0,T'),H''(R)nH-\R))nC\[0,T'),H''-^' 
with the same T' = T . See Yin fW/ for standard arguments. 
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By using the local well-posedness result in Lemma [T] and energy estimates, we obtain a 
precise blow-up scenario of the solutions to the Cauchy problem (I2.ip . 

Lemma 2. Let uq £ ff^(M) n H^^{M.) and u{t) be a solution of the Cauchy problem 
\2. 1\) in Lemma [IJ The solution blows up in a finite time T £ (0, oo) in the sense of 
liuitiT \\u{-,t)\\j{2 = oo if and only if 

lim sup u{x,t)ux{x,t) = +00. 



Proof. We only need to prove the necessary condition, since the singularity in u{x, t)ux{x, t) 
as i I T implies the singularity in t)||j:^2 as t | T. Assume a finite maximal existence 
time T £ (0, 00) and suppose, by the contradiction, that there is M > such that 

(2.3) supu{x, t)ux{x, t) < M < 00, VtG[0,T). 

Applying density arguments, we approximate the initial value uq G //^(M) by functions 
Uq £ n > 1, so that lim„^oo ""o = uq. Furthermore, write n"(t) for the solution of 

the Cauchy problem ()2.ip with initial data Uq . Using the regularity result proved in Lemma 
□ for s = 3, it follows from Sobolev's embedding that, if n"(t) £ C{[0,T),H^{R)nH-'^{R)), 
then u"'{x,t) is a twice continuously differentiable function of x on R for any t £ [0,T). It 
is then deduced from the short-pulse equation (jl.ip that 



di 



Kfdx = I u'^iulYdx < M / {u'iYdx 



and 

ulxfdx = 5 1 u'^UxiUxxfdx <5M I [u^^^dx. 



d_ 
lit 

The Gronwall inequality yields for all t £ [0, T), 



Since ||uq||j:/2 converges to ||uo||/f2 as n 00, we infer from the continuous dependence of 
the local solution u{t) on initial data uq that t)||j|^2 remains bounded on [0, T) for the 
solution u{t) in Lemma [TJ Therefore, the contradiction is obtained and either T is not a 
maximal existence time or the bound (j2.3p is not valid on [0,T). □ 

We also show that the blow-up of Lemma [2] is the wave breaking in the sense of condition 
(jl.3p . In other words, both ||ii(-, t)||/,oo and are uniformly bounded for all 

t£[0,T). 

Lemma 3. Let uq £ n H^^{R) and T > be the maximal existence time of the 

solution u{x,t) in LemmaUi Then, 

(2.4) M-,t)\\L^<Fi, \\d-\{-,t)\\L^ <Fo, t£[0,T), 
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where 
(2.5) 

(2.6) 



Fn 



l=yE, + E., + -E,Ft 



1/2 



Proof. For all t G [0,T) and the solution u{x,t), we have 



< 



X POO 

UUxdx — I 

-OO J X 

\u\\Ur\ 



UUxdx 



1 + yjl + uldx 



u% 



As a result, we obtain 



< eI'^ (^j u^{2+ ^Jl + ul- l)dx^ 

< i?i'/'(2i?o + i?ilK,t)||ioo)'/'. 

|n(-,t)||ioo <2So^i+^2||^(.^^)||2^^ 



1/2 



so that bound (|2.5p is found from the quadratic equation on ||u(-, i)|||oo . On the other hand, 
we have 

\\d^M;t)fm = \H-Mh + \\dxM-Mh 

= Eo + E_i + ^\\u{;t)tL^ 

< ^o + ^-i + ^^o||n(-,t)||ioo. 
Using the Sobolev inequality ||3~^ti||L=o < -^USar^nH/^i and bound (|2.5p . we obtain bound 

Let us introduce a continuous family of characteristics for solutions of the Cauchy problem 
(j2.ip . Let ^ G M, t G [0, T), where T is the maximal existence time in Lemma [U and denote 

x = X(e,t), u{x,t) = U{i,t), d-\{x,t) =G{^,t), 



(2.7) 
so that 

(2.8) 



X{t) = - 

x(o) = e 



1772 



f>(t) = G, 
C/(0) =txo(0, 



where dots denote derivatives with respect to time t on a particular characteristics x = 
X{^,t) for a fixed ^ G M. Applying classical results in the theory of ordinary differential 
equations, we obtain two useful results on the solutions of the initial- value problem ()2.8p . 
Conserved quantities E-i and Eq of the Cauchy problem (|2.ip are used to control values of 
U and G on the family of characteristics. 
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Lemma 4. Let uq G H'^(U) n H ^(M) and T > be the maximal existence time of the 
solution u{t) in Lemma d Then there exists a unique solution X{^,t) £ C^(R x [0,r)) 
to the initial-value problem i2.8\) . Moreover, the map X(-,t) : M M is an increasing 
diffeomorphism for any t G [0, T) with 

d^X{^, t) = exp u{X{i, s), s)u^{X{i, s), s)ds^ > 0, te[0,T), ^eR. 



Proof. Existence and uniqueness of X{^, t) G C^(]R, [0, T)) follows from the integral equation 

t) = ^-l f uHc, s)ds, t G [0, T), e e R, 

^ Jo 

since U{^,t) G C(M, [0,T)) for the solution u{t) in Lemma [TJ Using the chain rule, we 
obtain 

t) = -W{^, t)d^X{^, t) d^Xit t) = exp 1^ W{^, s)ds^ , 

where W{C,t) = u{X{i,t),t)u^{X{i,t),t) G C7(M, [0,T)). Therefore, d^X{i,t) > for all 
t G [0, T) and ^ G □ 

Let 

V{i,t) = u,{X{^,t),t), W{tt)=u{X{^,t),t)u^{X{^,t),t) ^ U{C,t)V{^,t) 
and compute their rate of changes along the family of characteristics 



(2.9) 



V = vw + u, 

W = + VG + U^. 



Let Fq, Fi > he fixed in terms of conserved quantities E^i, Eq, and Ei as in Lemma [3] 
and assume that there exists .fo ^ such that W{^q,0) > and 

/ 772 X 1/3 / N 1/2 

either 1^(^0,0)1 > ( , \V{Co,0)\W{^o,0) > Fi + (2Fo\V{^o,0)f - -Ffj , 

or 1^(^0,0)1 < (J^r^ ' , |l^(eo,0)|T^(eo,0) > Fi. 

Because of the invariance of system ()2.9p with respect to 

G ^ -G, U ^-U, V ^ -V, W ^W, 

it is sufficient to consider the case with V{^o, 0) > 0. We will prove that, under the conditions 
above, V{S,o,t) and W{S^Q,t) remain positive and monotonically increasing functions for all 
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t > 0, for which they are bounded, so that V{(,o, t) and W{£,Q,t) satisfy the apriori differential 
estimates 



(2.10) 



V > VW-Fi, 
W > W^-VFo. 



In what follows, we use V{t) and W{t) instead of V{(,o, t) and W{£,q^ t) for a particular ^ 
M. The following lemma establishes sufficient conditions on the initial point (y(0),T^(0)) 
that ensure that a lower solution satisfying 



(2.11) 



V = VW-Fi, 
W = W^-VFo, 



goes to infinity in a finite time. 

Lemma 5. Assume that the initial data for system 12.11\) satisfy 

72 \ 1/3 / T N 1/2 



(2.12) either V{0) > (^^) , ^(0)1^(0) > Fi + (^2FoVHo) - ^Ff^ , 

(2.13) or < V{0) < (^A.^ , V{0)W{0) > F^. 



Then the trajectory of system 112.11]) blows up in a finite time G (0, oo) such that V{t) 
and W{t) are positive and monotonically increasing for all t E [0,t*) and there is C > 
such that 

(2.14) limy(t) = oo, limW^(t) = oo, and lim{t^ - t)V{t) = C. 
Moreover, t^, is hounded by 

(2.15) t. < - ^(»' 



min { V{0), ( 1/2(0) - 2Foy3(0) + ^Ff ' ^''^ 



Proof. Let us first consider the homogeneous version of system (|2.11|) for Fi = 0, that is 
(2.16) 



V = VW, 



W = W^-VFn. 



Of course, Fi is never zero, otherwise Ei = 0. This case is used merely for illustration, 
since explicit solutions can be obtained for Fi = 0, whereas qualitative analysis has to be 
developed for Fi / 0. System (|2.16|) is integrable since 

where C = W{0)/V{0). Integrating the last equation for V{t), we obtain the explicit 
solution of the truncated system, 

Vit) = 7^^^, r^, W{t) = {C- Fot)V{t). 

l-CV{0)t + ^FoV{0)t^ ^ u; w 
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The solution reaches infinity in a finite time G (0, oo) if ^(0) > and C^V{0) > 2Fo. 
Note that these conditions coincide with condition (I2.12|) for Fi = 0. Also note that is 
the first positive root of 1 - CV{0)t + ^FoV{0)t'^ = so that U{Q = C- FqU > and 



^ _ W{0) - ^/W^jO) - 2FoV{0) 2^(0) ^ ^(0) 

V{0) + ^y2(0) -2Foy3(o) ~ ^Jv^{0) - 2FoV^O) 

Consider now the full system ()2.1ip . Let 

1 V 
V = -, W = - 

X X 

and rewrite the system in the form 

(2 17) / ^ = " ^' 

^ ' \y = Fixy-Fo. 

Expressing y from the system, we can rewrite it in the scalar form 
(2.18) y = Fix^ -X => x = f{x) + 3Fixx, 

where f{x) = Fq — F^x^. The only critical point of equation (j2.18p is (x, x) = [xq, 0) where 
xo = {Fo/FfY/^ is the root of f{x). 

We need to show that there is a domain C in the phase plane so that all 

initial data in D generate trajectories in D that cross the vertical line x = in a finite time. 
To do so, we construct a Lyapunov function for equation (j2.18p in the form 

E{x, x) = - Fox + ^Ffx^. 



The function E(x,x) has a global minimum at (xo,0). For any solution x{t), we have 

-E{x, x) = SFixx^ > for x > 0. 



dt 

The zero level of the Lyapunov function E{x, x) = passes through the points (0, 0) and 
(x=K,0), where x^, = A^^^xq > xq (see Figure [T]). It is clear that E{x,x) > in the domain 

D = {(x, x) : X > 0, X < o"(x)} , 

where 

cj(x) 



(2Fox- xG(0,x,), 
0, X G [x*, oo) 



We note that the condition (x, x) £ D is equivalent to the condition 

fo<x<x*, y> Fix2 + (2Fox - ^F^x^)^^^ , 
I X > X*, y > Fix^, 
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which is nothing but the set of conditions (|2.12p and (|2.13p at t = 0. By continuity, if 
(x, x) £ D at t = 0, then (x, x) remains in D for some time t > 0. 

No critical points of system ()2.18p are located in D and x{t) is decreasing function for 
any t > as long as the trajectory stays in D. Recall that E{x, x) > and ■^E{x, x) > 
for any £ D. A trajectory in D can not cross x = a{x) because E{x,a{x)) = for 

X £ (0, x*) and x = f{x) +?,Fixx < for x > x* and x < 0. Therefore, the trajectory either 
reaches x = in a finite time t^, G (0, oo) or escapes to x = — oo for x > 0. To eliminate the 
last possibility, we note that 

d f . 3Fi 2 



so that 



x{t) > x(0) - ^x^{0) + i/(x(0)) > -C30 

for any finite time interval. Moreover, x is bounded from zero in D by the level curve 
E(x,x) = E{x{0),x{0)), which is a convex curve in D. Therefore, 

x{t) < max{x(0),p}, t > 0, as long as x > 0, 

where p < is uniquely found from E{0,p) = E{x{0), x{0)), that is from the point of 
intersection of the level curve of E{x, x) = £'(x(0), x(0)) with the negative x-axis. Therefore, 

x{t) < x(0) +tmax{x(0),p}, 

so that x(t) reaches x = in a finite time t* G (0, cxd) for any trajectory in D. Moreover, 
finding p explicitly gives the bound on the blow-up time 

' n ^ ^ -f < ^^^^ 

U < X < X^, ^ (i2(o)_2Fox(0)+|F2x4(0))l/2 ' 

, ^ x{0) 



|i(0)|' 

which becomes bound (j2.15p after the return back to variable V{t). 

Since V = x~^, W = yx"^, and y = Fix"^ — x > 0, we have lim^it,^ = oo and 
lim^i-^^ W{t) = oo for any trajectories in D. Since x < for the trajectory in D, we also 
have x(t) ~ (t* — t) as t — > t=K so that there exists C > such that limf|f^(t* — t)V{t) = C. 
It remains to show that V{t) and W{t) are monotonically increasing functions on [0,t*). To 
do so, we write 

V = VW-Fi = -^, 
x^ 



x^ 



where 

g{x, x) = x^ — 2Fix^x — x/(x). 
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Figure 1: Domain D in the phase plane (x,x) of equation (j2.18p . 

For any trajectory in L>, x{t) < so that V{t) > 0. Furthermore, since g{x,x) is zero at a 
curve outside the domain D, because 

I 

g{x,a{x)) = Fox + -F^x"^ + 2Fix'^\a{x)\ > 0, xe {0,x^], 
then W{t) > for any trajectory in D. □ 



Recalhng that W = UV in system (j2.9p and using the bound ()2.4p . we obtain the upper 
bound for any solution at the family of characteristics 

(2.19) V = UV^ + U < FiV^ + Fi. 

We can now show that any upper solution with V^(0) > goes to infinity in a finite time. 
Lemma 6. Consider 

(2.20) V{t) = FiV^ + Fi 

with V{0) > 0. There exists G (0,oo) such that V{t) is positive, monotonically increasing 
for all t E [0, t^,) and there exists C > such that 



lim(t* - t)Vit) = C. 



Moreover, t^ < l/W{0). 



Proof. Since V > for any y € M, V{t) is monotonically increasing function. To show that 
V{t) reaches cxd in a finite time, one can integrate the separable equation ()2.20p explicitly 
and obtain 

V{t) = tan(arctan y(0) + Fit), 

so that 

_ 7r/2 - arctanF(O) ^ 1^1 

* ~ ¥i - FiV{o) - wxo)' 

since sup^jgi^ a;cot~"'^(x) < 1. □ 
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Applying results of Lemmas [5] and [H we conclude the proof of Theorem [3l 

Proof of Theorem 0. Let (V, W) satisfy system (j2.9p corresponding to the characteristics 
with ^0- Let (V, W) be the lower solution of system ()2.1ip in Lemma [5] with V_{0) = V{0) 
and W_{0) = W{0). Let V be the upper solution of equation (j2.20p in Lemma [6] with 
V{0) = V{0). Let be the blow-up time of the lower solution and be the blow-up time 
of the upper solution. 

The upper bound for the solution of system (|2.9p follows from the comparison principle 
for the differential equations since 

\vw + u\ = + i)\u\ < (y2 + i^p^ 

which implies that V{t) < V{t) for all t £ [0,t*), for which V{t) remains bounded. 
To obtain the lower bound, we note that 



V >V, W >W 
Let V = [V, Wf, V = [V,Wf, and 



V > VW -Fi> VW -Fi = V, 
W + FqV >W'^>W^ = W + FqV. 



L 




Fo 



be a nilpotent matrix of order one, so that e*^ = / + tL. Thus, we write 

V>V, W>W =^ 4- > 4- (e*-^V) . 

Integrating this equation in t > 0, we infer that 

e*^V(t) > e*^V(t), 
and since e*'^ is invertible for any t G M, we conclude that 

V{t)>V{t), w{t)>w{t), 

for all t G [0, r) C [0,t^), for which (y, VF) remain finite. Therefore, (F, VF) become infinite 
as t t T and T € □ 

Remark 3. The hounds ont^ andt^, in Lemmasl^ andl^i are inconclusive to compare T with 
the exact time of blow-up Tq := along a particular characteristic of the dispersionless 
advection equation 

3 Wave breaking in a periodic domain 

Consider now the Cauchy problem for the short-pulse equation (jl.ip in a periodic domain 

, { ut = ^u'^Ux + 5~^n, X G§, t > 0, 

\ u{x,0) = Uo{x), a; G S, 
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where S is a unit circle equipped with periodic boundary conditions and ^ is the mean-zero 
anti-derivative in the form 

px p px 

d~^u := / u{x',t)dx'— / / u{x' ,t)dx'dx. 
Jo Js Jo 

Local well-posedness and useful conserved quantities for the Cauchy problem (|3.1|) in a 
periodic domain are obtained in the following lemma. 

Lemma 7. Assume that uq S H^(S), s > 2 and J^uq(x) dx = 0. Then there exist a 
maximal time T > such that the Cauchy problem \3.1\) admits a unique solution 

u{t) G C{[0,T),H%§))nC\[0,T),H'-\S)) 

satisfying u{x, 0) = Uo{x) and /g u{x, t)dx = for all t € [0, T). Moreover, the solution u{t) 
depends continuously on the initial data uq and the quantities 

Eq = u^dx, El = y/l + uldx 
Js Js 

are constant for all t £ [0,T). 

Proof. Existence of the solution u{x, t) and continuous dependence on uq is proved on S 
similarly to what is done in Lemma [1] on M. To prove the zero-mass constraint, we note 

ut{t) G C{{^,T),W-\S)), u^uS) G C{%T),H'~\^)), 

so that for all t G (0, T), we have 

u{x, t) dx = i Utx dx -\ — / {u'^Ux)x dx = 0. 
Jn 2 j§ 

Initial values of Eq and Ei are bounded if uq G H'^{S), s > 2. Conservation of Eq and Ei 
on [0, T) follows from the balance equations 



X 



thanks to the continuity and the periodic boundary conditions for ^u{t) G C((0, T), if^+^(S)), 
u{t) G C{{0,T),H'{§)), and n^(t) G C{{0,T), H'-^{S)) in x on S if s > 2. □ 

Remark 4. The assumption J^Uo{x) dx = on the initial data uq in the periodic domain 
S is necessary as it follows from the following apriori estimate 

< - no||L2(s), VtG(0,r). 



u{x,t)dx— / uq{x) dx 
Js 

Note that f^u{x,t)dx = 0, for all t G (0,r) and u{t) G C([0, T), i?*(S)), s > 2. Hence the 
above estimate implies that J^uq{x) dx = 0. Note that no zero-mass constraint is necessary 
on an infinite line in Theorem{^ 
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The blow-up scenario for the solutions to the Cauchy problem (j3.ip coincides with the 
one in Lemma [2] after the change of M by S. The main result of this section is the proof of 
the finite-time wave breaking in a periodic domain, according to the following theorem. 

Theorem 4. Let uq £ H'^{§) and J^uo{x) dx = 0. Assume that there exists xq £ M. such 
that uo{xo)uq{xo) > and 



either 



or 



|no(xo)| > 
\uo{xo)\ < 



El 



1/3 



4K 



1/2 



, \Uo{Xq)\\U, 



f 1 \ 



El 



1/3 



4K 



1/2 



ho(2;o)lko(3;o)|^ > Ex. 



Then there exists a finite time T G (0, oo) such that the solution u{t) G C([0, T), i7^(S)) of 
the Cauchy problem ^3. 1\) blows up with the property 



limsupM(x,t)ti2^(x,t) = -|-oo, while \\m.\\u{-,t)\\L°° < Ei. 

Proof. Let T > be the maximal time of existence of the solution u(t) G C([0, T), i7^(S)) 
to the Cauchy problem (j3.ip constructed in Lemma [71 Since Jg u{x., t)dx = 0, for each 
t G [0, T) there is a .^^ G [0, 1] such that ^(^f , t) = 0. Then for x G S and t G [0, T), we have 



\u{x,t)\ 



Ux{x, t) dx 



it 



< 



/ \ux{x, t)\dx < El. 



Since d~^u{t) G C{[0,T), H^{S)) is the mean-zero periodic function of x for each t G [0, T), 
there exists another G [0, 1] such that d~^u{£,t,t) = 0. Then for x G S and t G [0,r), we 
have 



\dx^u{x,t)\= / u{x,t)dx < 

Jit 

Therefore, bounds (j2.4p are rewritten with 

Fq := ^/Eq, Fi := Ei 
The rest of the proof follows the proof of Theorem [31 



u{x, t)\dx < \/Eq. 



□ 



4 Numerical evidence of wave breaking 

The goal of this section is to complement the analytic results by several examples and 
numerical computations. More specifically, we first show that the sufficient condition for 
wave breaking in Theorem [3l is not satisfied for the exact modulated pulse solution to the 
short-pulse equation which is known to remain globally bounded in space and time. Then we 
consider the interplay between global well-posedness and wave breaking of Theorems [21 and 
[31 for a class of decaying data on an infinite line. Finally, we perform numerical simulations 
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in a periodic domain for a simple harmonic initial data and thus give illustrations to the 
sufficient condition for wave breaking in Theorem HI 

Theorem [3] gives a sufficient condition for formation of shocks in the short-pulse equation 
(jl.ip on the infinite line. Let us show that this condition is not satisfied for exact modulated 
pulse solutions obtained in [U [lO] . The simplest one-pulse solution is given in the parametric 
form 

u{x,t) = U{y,t), x = X{y,t), 



where 



(4.1) 



U{y, t) = Amn 



X{y, t) = y + 2mn 



m sin ijj sinh (j) + n cos ip cosh i 
m? siv? ip -\- n'^ cosh^ (j) 

m sin 2'tp — n sinh 2(j) 



sin^ ^lJ + n"^ cosh^ cj) ' 
m G (0, 1) is an arbitrary parameter, n = Vl — Tn"^, and 

(p = m{y + t), V = ^{y ~ i)- 
The pulse solution enjoys the periodicity property 
r U{y,t) = U{y-^,t + ^), 
I X{y,t)=X{y-^,t + ^)+^ 
and an exponential decay in any direction transverse to the anti-diagonal on the {y, t)-plane. 
Since 



dX 
dy 



Sm'^n? sin'^ -i/^ cosh'^ 



1 



(m^ sin^ i/j + n'^ cosh^ (f))'^ 
the function x = X{y, t) is invertible in y for all t G M if 

msin^ 



cos I 4arctan 



msinip 
n cosh 6 



n cosh < 



< tan ■ 



m vr 
— < tan — . 

n - 8' 



that is for all m £ {0,mcr), where rricr = sin ^ ~ 0.383. For these values of m, the pulse 
solution u{x, t) is analytic in variables (x, t) G M, has the space-time periodicity 



u{x, t) 



u X 



IT IT 
-,t + - 

m m 



{x,t) G 



and the exponential decay in the transverse direction to the anti-diagonal in the (x, f)-plane. 
The graph of a nonsingular pulse solution for m = 0.32 is shown on Figure [2] (left). 



Remark 5. Coordinate y in the exact solution j[ ) is different from coordinate ^ in the 
method of characteristics because X(y,0) ^ y. Nevertheless, X{y,t) and U{y,t) satisfy the 
same set of equations 

dX 1 , dU 

so that ^ and y are uniquely related by the representation ^ = X{y,0). If y is found as a 
function of E,, the initial data of the Cauchy problem ^2.1\) is found from tto(0 = U{y,0). 
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Since u{x, t) is analytic in a; G M for any fixed t G M and decays to zero exponentially fast 
at infinity, it is clear that u{-,t) G H'^{M.). Furthermore, since d~^u = ut — ^u^Ux, it also 
follows that G H~^{M.). We compute numerically bounds Fq and Fi using the exact 

solution ()4.ip and new definitions 



sup||u(-,t)||j 



Fo := sup \\dx ^u{-,t) 



It follows from Remark [T] that bounds Fi and Fq defined above preserve the sufficient 
condition of Theorem [3] with respect to the scaling transformation (jl.Sp . 



Let us define 



/i := sup [|no(x)||uo(x)p - Fi] , 



(4.2) 



where 



/2 := sup 



1 \ 1/2' 

\uo{x)\\u'o{x)\' - Fi - ( 2Fo|n'o(x)|3 - -Ff 



X G M : \u'o{x)\ < (^^^ , uo(x)n'o(x) > 0, | , 



X G 



/ Ff \ 



\u'oix)\ > [j^j > 'Uo(a;)no(x) > 0, 



According to Theorem [3l wave breaking occurs if either fi or /2 is positive. For the exact 
modulated pulse solution (j4.ip at t = 0, the numerical calculations show that the set I2 is 
empty and the quantity fi is strictly negative for any m G (0,mcr), see Figure [2] (right). 
Therefore, the sufficient condition for the wave breaking in Theorem [3] is not satisfied, 
which corresponds to our understanding that the exact modulated pulse solutions (|4.ip 
remain bounded for all {x,t) G M?. We note, however, that the sufficient condition for 
the global well-posedness in Theorem [2] is satisfied only for pulses with m G (O, ■^), since 
2^/2ElE2 = 32m. This computation shows that the sufficient condition of Theorem [2] is 
not sharp. 

Next we compare the sufficient conditions for the global well-posedness and wave breaking 
in Theorems [2] and [3] for a class of initial data 



(4.3) 



uo{x) = a{l 



2bx^)e 



a > 0, 6 > 0, 



where parameters a and b determine the amplitude and steepness of uq. Note that the 
zero-mass constraint (j2.2p is satisfied by uq and it is clear that uq G -fr^(M) H H^^(R). The 
conserved quantities F_i and Eq can be computed analytically, 



F_i 



'^7^(256^/2 - 51a26) 



2048 V63 



En 



3a2\/27f 



whereas the conserved quantities Ei and F2 are not expressed in a closed form. Using 
numerical approximations of the integrals, we determine the boundary of the well-posedness 
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Figure 2: The exact modulated pulse solution (j4.ip of the short-pulse equation (jl.ip for 
m = 0.32 (left). The quantity /i is negative for any m E (0, mc^) (right). 

region in the (a, 6)-plane by finding the parameters a and b from the condition 2^2EiE2 = 1. 
We also compute the boundary of the wave breaking region in the (o, 6)-plane by computing 
fi and /2 in (j4.2p . Unlike the case of modulated pulses, we find that the set Ji is empty 
and /2 may change the sign along the curve on the (a, 6)-plane. The two boundaries are 
shown on Figure [31 where we can see that the two regions of global well-posedness and wave 
breaking are disjoint. 

Finally, we perform numerical simulations of the periodic Cauchy problem ()3.ip with the 
1-periodic initial data 

(4.4) lio(x) = a cos 27rx, a > 0. 

The two conserved quantities Eq and Ei in Lemma [7] are computed analytically as 

Eo = ia2. El = -E(27rai), 
z vr 

where E stands for a complete elliptic integral. Using the above conserved quantities we 
find out that the sufficient condition for the wave breaking in Theorem S] is satisfied for 
a > 1.053. 

To illustrate the behaviour of a solution u{x, t) to the Cauchy problem (j3.ip . we perform 
numerical simulations using a pseudospectral method. When the parameter a is sufficiently 
small, the value of 

W{t) := sup u(x,t)ux{x,t) 

remains bounded as shown on the top panel of Figure H] for a = 0.2. On the other hand, 
when a becomes larger, the wave breaking occurs, as on the bottom panel of Figure H] for 
o = 0.5. On the bottom right panel of Figure H] we show using the linear regression that 
the curve W~^{t) is fitted well with the straight line A + Bt for some coefficients {A,B). 
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0.5 1 1.5 2 



Figure 3: Boundaries of the global well-posedness and the wave breaking in the Cauchy 
problem (12. 1[) with initial data (|4.3p : the global well-posedness occurs below the lower 
curve and the wave breaking occurs above the upper curve. 



0.18 





Figure 4: Solution surface u{x,t) (left) and the supremum norm W{t) (right) for a = 0.2 
(top) and a = 0.5 (bottom). The dashed curve on the bottom right picture shows the linear 
regression with C = 1.072, T = 1.356. 
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Figure 5: Time of wave breaking T versus a (left). Constant C of the linear regression 
versus a (right). 

Thus, we make a conclusion that 

W{t) ~ - — - for < T - t < 1, 

where C = —B^^ and T = —AB^^. Using the linear regression, we also obtain pairs (T, C) 
for different values of a. The results are shown on Figure [5l Note that the constant C 
approach 1 as a gets larger. This observation is consistent with the exact blow-up law 
= obtained for the dispersionless advection equation (jl.2p using the method of 
characteristics. 
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